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ABSTRACT 

The smaller the angular scales on which the anisotropies of the cosmic microwave background 
(CMB) are probed the more important their distortion due to gravitational lensing becomes. 
Here we investigate the maxima and minima of the CMB lensing deflection field using gen¬ 
eral extreme value statistics. Since general extreme value statistics applies to uncoii'elated 
data in first place we consider appropriately low-pass filtered deflection maps. Besides the 
suppression of correlations filtering is required for another reason: The lensing field itself is 
not directly observable but needs to be (statistically) reconstructed from the lensed CMB by 
means of a quadratic estimator. This reconstruction, though, is noise dominated and therefore 
requires smoothing, too. In idealized Gaussian realizations as well as in realistically recon¬ 
structed data we find that both maxima and minima of the deflection angle components follow 
consistently a general extreme value distribution of Weibull-type. However, its shape, location 
and scale parameters vary significantly between different realizations. The statistics’ potential 
power to constrain cosmological models appears therefore rather limited. 
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1 INTRODUCTION 

While traversing the Universe photons of the Cosmic Microwave 
Background (CMB) are deflected hy the intervening large scale 
structure: They are gravitationally lensed. Lensing widens the size 
distrihution of hot and cold spots imprinted at the surface of last 
scattering. At the primary acoustic peaks this effect amounts to sev¬ 
eral per cent and dominates in the damping tail transferring power 
from larger to smaller angular scales. CMB lensing has been mea¬ 
sured in a number of experiments, most recently in the Planck data 
(Planck Collaboration 2014b). 

Precise knowledge of the lensed CMB temperature and po¬ 
larization power spectra allows to lift several parameter degenera¬ 
cies (eg. Smith et al. 2006). In particular the projection degeneracy 
of the primary CMB between spatial curvature and a cosmolog¬ 
ical constant can be alleviated through lensing (Hu et al. 2006). 
Moreover, CMB lensing is expected to contain important informa¬ 
tion about the onset of late-time cosmic acceleration (Acquaviva 
& Baccigalupi 2006) and has recently been reconsidered in the in¬ 
vestigation of early dark energy (Pettorino et al. 2013). Regard¬ 
ing the statistics of the observed, i.e. lensed, CMB anisotropies 
there are two important differences with respect to the primordial 
fluctuations. Firstly, the lensed CMB is necessarily non-Gaussian, 
even in case of Gaussian temperature fluctuations, because lensing 
is a nonlinear effect, and secondly it is not statistically isotropic. 
The off-diagonal elements of its covariance can be used to recon¬ 
struct (in a statistical sense) the lensing deflection field (Hu 2001; 
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Okamoto & Hu 2003; Hirata & Seljak 2003) and the correspond¬ 
ing power spectrum (Kesden et al. 2003; Hanson et al. 2011). The 
non-Gaussian structure of the observed CMB due to lensing has 
been directly addressed by Hu (2000) and Lewis et al. (2011) fo¬ 
cusing on the CMB bispectrum, while Munshi et al. (2013) em¬ 
ployed skew-spectra and Minkowski Functionals. The reconstruc¬ 
tion of the CMB deflection field is of great importance in removing 
lensing induced B-mode polarization in order to access primordial 
tensor perturbations (Knox & Song 2002; Seljak & Hirata 2004). 
Furthermore, since the deflection field probes the integrated mass 
distribution between today’s observer and the surface of last scatter¬ 
ing it helps inventory the Universe’s matter content (Lesgourgues 
et al. 2006). 

All these applications of CMB lensing depend crucially on an 
accurate theoretical modeling of the lensing effect. Most work in 
this field concentrated on the exact computation of lensed power 
spectra (Seljak 1996), culminating in the all-sky expressions of 
Challinor & Lewis (2005). In addition to that, there are investi¬ 
gations of the corresponding two-point correlation function (e.g. 
Metcalf & Silk 1997; Takada et al. 2000). The deflection field it¬ 
self received comparatively little attention. Carbone et al. (2008) 
constructed realistic CMB deflection maps from the Millennium 
Simulation. The resulting lensed CMB temperature and polariza¬ 
tion maps agree very well with the theoretical prediction (Carbone 
et al. 2009). Using the results of Carbone et al. (2008), Merkel & 
Schafer (2011) showed that the small amount of non-Gaussianity 
present in the CMB lensing excursion angle has only a marginal 
effect on the lensed CMB spectra. In a recent study Carbone et al. 
(2013) produced lensing deflection maps from A-body simulations 
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in the framework of interacting dark energy cosmologies. They 
found that for coupled dark energy models the (statistical) power of 
lensing distortions are either enhanced or suppressed (depending on 
the specific model under consideration) at the ten per cent level with 
respect to a A cold dark matter (ACDM) universe. Such differences 
are not necessarily restricted to the two-point function but may al¬ 
ready be present at lower level. A particularly powerful one-point 
statistics is provided by the distribution of extreme values. Distinct 
properties of such a distribution then might be used to discriminate 
between different cosmologies. Furthermore, extreme value statis¬ 
tics are potentially more sensitive to deviations from Gaussianity 
than conventional higher-order statistics because the distribution of 
peaks in a Gaussian field declines rather fast. Non-Gaussianities, 
however, alleviate this sharp fall-off. In this paper we investigate 
to which extent and under which conditions general extreme value 
statistics can be applied to the CMB deflection field in the context 
of standard ACDM. 

The outline of this paper is as follows: In Section 2 we intro¬ 
duce the basic concepts of general extreme value statistics and ex¬ 
plain its application to CMB lensing. Section 3 compiles our results 
for both simulated Gaussian and reconstructed data. We conclude 
in Section 4. 

As reference we choose a ACDM cosmology with parame¬ 
ter values: Ocdm = 0.1143, = 0.02256 and Da = 0.721 

describing the Universe’s energy budget; = 2.457 x 10“^ and 
ris = 0.96 characterizing the primordial perturbations; h = 0.701 
for today’s Hubble constant, i.e. Hq = 100 km s“*. 


2 METHODOLOGY 

2.1 CMB lensing 

Due to Liouville’s theorem gravitational lensing conserves the sur¬ 
face brightness of a given source (e.g. Bartelmann & Schneider 
2001). Thus, the lensed temperature contrast 0 results from the pri¬ 
mordial one via the remapping 

@(«) = Q[h + </(«)] ( 1 ) 


(see Lewis & Challinor 2006, for a comprehensive review on CMB 
lensing). The lensing deflection angle d is the angular gradient of 
the lensing potential 


d(n) = V 0 , 0 i/r(n) = 2 



XX* 


( 2 ) 


where the latter follows from the line-of-sight projection of the 
gravitational potential T*. Assuming a spatially flat universe we can 
use angular diameter and comoving distances interchangeably (de¬ 
noted hy x)- We further take recombination to be instantaneous, i.e 
the last scattering surface at^* is infinitely thin. As detailed in the 
appendix of Lewis (2005) the remapping of the primordial CMB 
on the full sky involves the two components of the deflection an¬ 
gle and its modulus d as well. In what follows we thus consider 
four different fields, both components of the spin -1 field d and the 
scalar fields if/ and d. 


2.2 CMB lensing reconstruction 

The nonlinear nature of CMB lensing becomes apparent from equa¬ 
tion (1). The statistics of the lensed temperature fluctuations is 
therefore inevitably non-Gaussian and its covariance matrix ac¬ 
quires off-diagonal elements which are, to lowest order, linear in 



Figure 1. Reconstruction noise of the quadratic estimator of the CMB lens¬ 
ing potential assuming a Planck-liks experiment (cf. Hu & Okamoto 2002). 
The noise levels correspond to different combinations of the CMB temper¬ 
ature (0) and polarization (£/B-modes). The ©©-estimator is close to the 
minimum variance combination (mv) of all estimators due to the limited po¬ 
larization sensitivity of Planck. The black dashed cuive indicates the power 
spectrum of the deflection field. 


the lensing potential. Thus, suitably weighted they can be used to 
construct a quadratic estimator of the lensing potential (Okamoto & 
Hu 2003). This estimator is unbiased and optimal in the sense that 
its (Gaussian) variance is minimal. However, as Figure 1 demon¬ 
strates for CMB experiments with sensitivity comparable to the 
Planck satellite the reconstruction is noise dominated. Furthermore, 
we see that the estimator based on the temperature field alone per¬ 
forms best. The gain from including polarization data is negligible. 
We therefore focus on the ©©-estimator and leave other combina¬ 
tions aside. Once the lensing potential is reconstructed the compo¬ 
nents of the deflection angle are readily obtained by differentiation. 


2.3 General extreme value distribution 


Under rather general conditions the extrema of a set of random vari¬ 
ates follow a distinct family of distribution functions (e.g. Gum- 
bel 1958; Castillo et al. 2004). To be specific, consider the max¬ 
imum M„ = max(Yi ,... ,X„) of a sequence of independent and 
identically distributed (i.i.d.) random variables A,. An analog to the 
central limit theorem states that in the limit n oo the cumulative 
distribution function (CDF) of M„ approaches the general extreme 
value distribution (GEV) given by 


Gaj3.y(x) = exp 



i/r 


(3) 


with location, scale and shape parameters a, jS and y, respectively. 
The sign of the shape parameter determines the type of the GEV: 
Gumbel-type (y = 0), Frechet-type (y > 0) and Weibull-type 
(y < 0). Analogous considerations hold true for the CDF of min¬ 
ima. The corresponding probability density function (PDF) results 
from differentiating equation (3) 


Saji.yix) — 


dGfl.^ y(x) 

dx 


(4) 


There are numerous applications of extreme value statistics in cos¬ 
mology (e.g. Waizmann et al. 2012; Capranico et al. 2013, consid¬ 
ering strong and weak lensing, respectively). A CMB related study 
has been carried out by Mikelsons et al. (2009), while Colombi 
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Figure 2. Illustration of our sampling strategy. A map of the CMB lensing 
potential ist covered by 1000 patches. The radius of each patch is 2°. 

et al. (2011) presented a general investigation of smooth Gaussian 
random fields in two and three dimensions. 

2.4 Sampling strategy 

In case of a collection of normal Gaussian variates the distributions 
of maxima and minima each approach the standard Gumbel-type 
CDF, i.e. Go, 1 , 0 . Both assumptions, Gaussianity and statistical inde¬ 
pendence, hold neither for the lensing potential nor the correspond¬ 
ing deflection field. The deviations from Gaussianity are small but 
the deflection field is correlated over the degree scale correspond¬ 
ing to the angular size of a typical lens which is located half the way 
to the last scattering surface (cf. Lewis & Challinor 2006). These 
correlations can be partially removed by applying a low-pass filter 
to the harmonic representation of the lensing potential 

= J'dQYfJfi)i/f(n). (5) 

Dividing the smoothed maps into patches we can estimate the PDF 
of extreme values by determining the maximum, respectively min¬ 
imum, in each patch. A faithful estimate put two major constraints 
on the choice of patches. On the one hand they have to be extended 
enough to contain a sufficiently large number of deflection angles 
in order to make the sample converge to the limiting CDF. On the 
other hand the patches must not overlap but have to be as numerous 
as possible to sample the PDF appropriately. In order to satisfy both 
requirements, we cover the smoothed map homogeneously with 
1000 discs the radius of which we choose to be 2°. This corresponds 
to a ~30 per cent coverage of the whole sky. The correlation coef¬ 
ficient between two patches has been estimated numerically from 
a sample of reconstructions and is typically ~0.3. An exemplary 
collection of patches is shown in Figure 2. The coordinates of the 
discs’ centers are drawn from a random distribution being uniform 

over §2. 

3 RESULTS 
3.1 Gaussian data 

We start our analysis with a sample of 100 Gaussian realizations of 
the lensing potential and the corresponding deflection field. For the 
generation of the synthetic data we first compute the lensing po¬ 
tential power spectrum using cams* (Lewis et al. 2000). Then 

* http://camb.info/ 


we invoke the healpix^ package (Gorski et al. 2005) to generate 
a set of harmonic coefficients such that = C'^'^6tc>6„„„>. 

Subsequently, these coefficients are used to build full-sky maps of 
the lensing potential and its first (angular) derivatives with healpix 
resolution parameter Aside = 2048, which captures angular scales 
of ~ 1.72'. During the map making procedure we include multipoles 
up to CS. = 3000. 

In the analysis we choose a sharp cut-off in the harmonic do¬ 
main (^max = 100) as filter function. Following the sampling pro¬ 
cedure outlined above we estimate the distribution function of the 
extrema for each realization. We fit the GEV (3) to the data al¬ 
lowing all three parameters to vary. Shape, location and scale pa¬ 
rameters are obtained by maximizing a log-likelihood function. A 
typical fit result is illustrated in Figure 8 for the maxima of the 
6 -component of the lensing deflection angle. Note that data (grey 
boxes) and fit (red dashed curve) in this figure are derived from a 
realistic reconstruction (which we discuss at length in Section 3.2) 
instead of a Gaussian realization. Maxima are generally denoted by 
a superscripted ’-l-’, e.g. for the lensing potential, while we use 
a superscripted ’-’in case of minima. 

We investigate the fit results in a more quantitative way in Fig¬ 
ure 3. Here we concentrate on the shape parameter. We find that for 
all four observables the distribution of maxima is consistently de¬ 
scribed by the Weibull-type GEV. The (relative) scatter, however, 
is large. It is largest in case of the lensing potential and the deflec¬ 
tion angle modulus, where it amounts to about 50 per cent. This 
compares to approximately 30 per cent in case of the deflection 
angle components. As expected from statistical isotropy the shape 
parameter is very much the same for the two deflection angle com¬ 
ponents (y ~ -0.24) taking the fit uncertainty into account. The 
shape parameter of the modulus, which is bound by zero from be¬ 
low, is markedly larger, i.e. less negative (y ~ -0.12), whereas we 
find a smaller value (y -0.29) for the lensing potential. Recover¬ 

ing a GEV of Weibull form agrees with the general considerations 
of Colombi et al. (201 1). 

Concerning the distribution of minima we find comparable re- 
sulfs. They are well described by the Weibull-type (see Figure 3) 
and the scatter is of the same order. As expected, the minima of the 
modulus do not follow the GEV due to its lower bound. The poor 
fit results of the shape parameter are omitted. 

In a next step we consider the influence of the cut-off choice. 
We reexamine the GEV of our sample for different filter scales. We 
aim at rather general properties which are independent of a particu¬ 
lar realization. Recalling the spread in shape parameters we believe 
that the deflection field components are most likely to comply with 
the requirement of generality. We therefore leave the other two ob¬ 
servables aside in the following. In addition to the degree-size cut¬ 
off we investigate the range = 500,..., 3000 in steps of 500. 
For each cutoff we plot the mean over the entire sample in Figure 4. 
The corresponding distributions functions are shown in Figure 5. 
The error bars in Figure 4 indicate the standard deviation which is 
at the 5-15 per cent level. Smallest (largest) scatter is observed in 
the location (shape) parameter (note that for the location param¬ 
eter the errors are artificially amplified by a factor of ten in the 
plot). Figure 4 reveals that the results for the shape parameter are 
rather insensitive to the cut-off choice. The other parameters of the 
GEV depend much more strongly on the filter scale. In particular 
we observe that removing small-scale fluctuations results in smaller 
location and scale parameters. Thus, mean and width of the corre- 
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Figure 3. Fit results for the shape parameter y of the GEV considering the 
maxima (’+’) and minima of the lensing potential (i/'*), the modu¬ 
lus (d*) and components (dg^ of the corresponding deflection field. The 
minima of the modulus do not follow a GEV function. The en'or bars cor¬ 
respond to the estimated Icr error of the fit, while the grey dashed line indi¬ 
cates the mean value. We only show half of the realizations for clarity. 


spending GEV decrease (cf. Figure 5). This behaviour is expected 
because the degree-size cut-off has been introduced to weaken the 
correlations in the lensing field. Its application removes inevitably 
the contributions of the rather efficient degree-sized lenses. Finally, 
we would like to emphasize that our findings presented so far do not 
show any significant differences for both our reference model and 
the Planck best fitting cosmology (Planck Collaboration 2014a). 
The main difference between these two parameter sets, which is 
particularly relevant for lensing, is that in case of Planck the mat¬ 
ter density is enhanced, whereas the amplitude of the primordial 
fluctuations is reduced. Considering the minima leads to the same 
conclusions. 

We next demonstrate that our results are independent of the 
patch size chosen for the sampling process explained in Section 2.4. 
In Figure 6 we show the fitting results for the shape parameter as 
a function of patch size. As before we use 100 maps of the deflec¬ 
tion field and lensing potential, respectively, which are smoothed 
in harmonic space at = 100. The disc radii range from 2° to 
8°, while the percentage of coverage is kept fixed to be 30 per cent. 
Thus, the larger the radius the smaller the number of discs covering 
the map in order to avoid overlapping patches. From Figure 6 it can 
be seen that the shape parameter hardly varies once the sampling 
process is able to resolve the limiting distribution. 


3.2 Reconstructed data 

Having found that the distribution of extrema in Gaussian realiza¬ 
tions of the CMB lensing deflection field is well fitted by a Weibull 
GEV we now leave these idealized conditions behind and address 
the more realistic case of reconstructed deflection fields. To this end 



f 
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Figure 4. Location a, scale and shape parameter y (from top to bottom) 
for the distribution of maxima of the deflection angle components as a func¬ 
tion of the cut-off in the harmonic domain. The error bars correspond to the 
Icr error over 100 Gaussian realizations and have been artificially amplified 
by a factor of ten in case of the location parameter. Filled symbols indi¬ 
cate our reference model, while light ones denote the Planck best-fitting 
cosmology. 
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Figure 5. Distribution functions of deflection angle maxima for different 
cut-off choices. The shifted curve for = 3000 has the same location 
parameter as that found for tmux = 100 and illustrates the broadening of the 
distribution with increasing cut-off size. Dashed lines indicate the Planck 
best-fitting cosmology. 


we generate 100 lensed CMB maps using lenspix'^ Lewis (2005) 
and subsequently apply the optimal quadratic estimator described 
in Section 2.2 to recover the lensing potential. In the framework 
of LENSPix lensed maps are generated by remapping the pixels of 
unlensed CMB maps in combination with a bicubic interpolation 
scheme. We choose the resolution of the unlensed maps to be a 
factor of 1.5 higher than that of the lensed maps for which we use 
HEALPix resolution Wide = 2048 and set = 3000. The initial 
CMB power spectmm from which we synthesize the input maps is 
computed using again cams. Note that lenspix uses a purely Gaus¬ 
sian realization of the lensing potential to construct the lensed CMB 
maps because it only considers the power spectrum of the poten¬ 
tial; no higher statistics are included. Hence, the lensing potential 


^ http://cosmologist.in£o/lenspix/ 
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disc radius [degree] 


Figure 6. Fit results for the GEV shape parameter y of the maxima of the 
CMB lensing potential and both deflection field components as a function of 
disc radius, i.e. patch size, used for sampling. En'or bars indicate Icr en'ors. 
Note the artificial olfset in radius we included for clarity in case of the 
deflection field components. 
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Figure 7. Fit results for the location, scale and shape parameters of the GEV 
considering the maxima of 100 reconstructed deflection angle maps. The 
error bars indicate an estimate of the fit uncertainty. Solid lines represent 
the average over all realizations, while the dashed ones correspond to the 
mean values found from Gaussian realizations. 


realizations considered here differ from those of Section 3.1 by the 
noise due to the estimator (cf. Figure 1). 

This noise makes a low-pass filter increasingly important. 
Since the deflection field is derived by differentiation from the 
lensing potential, there is a simple algebraic relation between 
the harmonic coefficients of the potential and those of the spin- 
1 components of the deflection field (e.g. Lewis 2005): idt„, - 

+ 1) if/cm- Thus, including fhe quite noisy small scale modes 
in fhe potenlial reconstruction introduces a substantial bias because 
the noise leads to a spurious enhancement of the deflections and 
thus widens the GEV. This enhances artificially both location and 
scale parameter. They increase by a factor of about 1.5-2. The 
shape parameter is only marginally affected, the distribution of 
maxima follows a Weibull-type GEV. We illustrate these findings 
in Figures 7 and 8 (left panel). The scatter observed in Figure 7 is 
compatible to the case of Gaussian realizations (cf. Figure 3). The 
right panel of Figure 8 reveals that in case of the reconstructed lens- 
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p-value of KS test 


Figure 9. Cumulative frequency of the outcome of Kolmogorov-Smimov 
(KS) test results. The distribution of maxima in the reconstructed CMB 
deflection angle components has been tested against a Weibull GEV the pa¬ 
rameters of which have been inferred from a different data set (thick curves). 
In addition we show the outcome for both maxima and minima in the re¬ 
constructed lensing potential (light curves). The grey boxes indicate the out¬ 
come for samples drawn from the assumed GEV. 


ing potential the discrepancy with respect to the mock data is even 
larger. The GEV of (both) potential maxima and minima obtained 
from Gaussian simulations is (are) much narrower than that(those) 
for the reconstructed data. 

These findings make, to say the least, a unified description 
of the maxima distribution in terms of the averaged shape, loca¬ 
tion and scale parameters rather questionable. In order to address 
our concerns quantitatively we investigate the goodness of fit em¬ 
ploying a Kolmogorov-Smirnov (KS) test. It is a non-parametric 
test and in order to get meaningful results the distribution param¬ 
eters must not be estimated from the same data set the test is ap¬ 
plied to. We therefore generate a second sample of 100 reconstruc¬ 
tions which are completely independent of the one presented before 
and from which the distribution parameters were determined. Both 
samples have only the background cosmology in common. The out¬ 
come of the test, i.e. the KS supremum values between proposed 
and fitted CDF, can be translated into a more readily accessible p- 
value. Usually, a p-value less than five per cent is regarded as an 
indicator that there is strong presumption against the null hypoth¬ 
esis, that the data follows the proposed distribution (in our case a 
Weibull GEV parametrized by the average shape, location and scale 
parameters), whereas an outcome greater than ten per cent implies 
no presumption against the null hypothesis. We compile the test re¬ 
sults in Eigure 9. Eor comparison we also show the results for 100 
data sets which have been sampled from the assumed Weibull dis¬ 
tribution. Considering the 6-component only ~ 40 (~ 25) per cent 
of the reconshucted maps exceed the threshold of 0.05 (0.1). The 
agreement is better for the other component, where ~ 75 (~ 55) per 
cent of the fit results pass the first (second) significance level. How¬ 
ever, this needs to be compared to the outcome of the mock data. 
Not more than ten per cent of it fall below a significance level of 
0.1. Thus, the distribution of maxima of the CMB deflection angle 
components is consistently described by a Weibull-type GEV. But 
we find no evidence for a universal set of shape, location and scale 
parameters independent of the particular realization. 

In order to investigate whether the lensing potential suffers 
from the same lack of universality we repeat the KS test for the 
extrema of the reconstructed potential. The results are shown in 
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Figure 8. Illustration of the PDF of maxima of the 6-component of the reconstructed CMB deflection angle (left panel) and the corresponding lensing potential 
(right panel). The data is well fitted by a Weibull GEV when its parameters are infen'ed from the data itself (dashed curve). The disagreement is apparent when 
the averaged shape, location and scale parameters obtained from either 100 Gaussian realizations (dashed dotted curve) or 100 reconstructions (solid curve) 
are used in stead. 


Figure 9, too. Only a tiny fraction of fits passes the test; roughly 
20 (10) per cent of the fit results are above the first (second) signifi¬ 
cance level. This is in line with our previous result that the (relative) 
scatter in the shape parameter is much larger for the lensing poten¬ 
tial than for the deflection field (cf. Section 3.1). 

Despite the missing universality we would finally like to 
roughly estimate under which conditions the GEV of CMB lens¬ 
ing angles might help find deviations from a ACDM cosmology. 
A shape parameter which exceeds clearly the typical fluctuations 

seen in Figure 7, say y -0.1 or -0.4 to be specific, would hint at 

a different cosmological model. Such a model, however, must have 
rather distinct lensing properties (compared to ACDM): Assuming 
that both location and scale parameters are unchanged the expec¬ 
tation value of the maximum deflection angle would become ten 
per cent larger or smaller, respectively. The variation is less pro¬ 
nounced when the scale parameter, too, changes at a comparable 
level. For the detection of such strong variations in the lensing ef¬ 
ficiency, however, cosmic shear measurements are much more ap¬ 
propriate. 


4 CONCLUSION 

We have applied general extreme value statistics to CMB lensing. 
We have found that both maxima and minima of the lensing poten¬ 
tial and the deflection angle components derived from it follow a 
GEV of Weibull form. This is also trae for the maxima of the de¬ 
flection field modulus. A basic assumption of general extreme value 
statistics is statistical independence. We therefore employed a sharp 
cut-off in harmonic space in oder to suppress correlations in the 
lensing potential. Low-pass filtering becomes particularly impor¬ 
tant when a realistic but noisy reconstruction scheme for the lens¬ 
ing potential, which is not directly observable, is employed. Since 
shape, location and scale parameters of the GEV vary considerably 
between individual realizations of the deflection field the statistics’ 
practical use as an additional cosmological probe is rather limited. 
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